Kinetics of shape equilibration for two dimensional islands 
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We study the relaxation to equilibrium of two dimensional 
islands containing up to 20000 atoms by Kinetic Monte Carlo 
simulations. We find that the commonly assumed relaxation 
mechanism - curvature-driven relaxation via atom diffusion 
- cannot explain the results obtained at low temperatures, 
where the island edges consist in large facets. Specifically, 
our simulations show that the exponent characterizing the 
dependence of the equilibration time on the island size is dif- 
ferent at high and low temperatures, in contradiction with the 
above cited assumptions. Instead, we propose that - at low 
temperatures - the relaxation is limited by the nucleation of 
new atomic rows on the large facets : this allows us to explain 
both the activation energy and the island size dependence of 
the equilibration time. 



I. INTRODUCTION 

There is a continued interest in the understanding, 
description and control of structures at the nanometer 
scales 10-0] . This is partially due to technological applica- 
tions of nanostructures, and partially to the fundamental 
interest of understanding how macroscopic concepts can 
(or cannot) be extrapolated down to these scales. On the 
technological side, a controlled preparation and conserva- 
tion of these structures demands a comprehension of their 
time evolution, which could be rapid due to the small 
scales involved (typical structures contain some hundred 
atoms). From the fundamental point of view, it is in- 
teresting to investigate how the theoretical tools which 
have been developed to deal with the kinetic evolution 
of macroscopic objects (size larger than a micrometer) 
by Herring, MuUins and Nichols |^, which are based on 
coarse-grained, continuous equation, can be used at the 
nanometer scale. For example, one could wonder whether 
sintering of ceramic or metallic nanopowders can be an- 
alyzed with these classic tools since it is not clear that 
macroscopic concepts such as curvature, chemical poten- 
tial, etc. should retain their relevance when dealing with 
structures containing only few atoms. 

Here, we focus in the dynamics of equilibration of two 
dimensional (2D) nanocrystallites starting in an out of 
equilibrium shape. MuUins' theory of shape relaxation 



is based on the curvature being well defined. Then, two 
questions arise. Can one use partial differential equations 
to study the shape relaxation of faceted nanocrystals? 
Can one use them at all far from the thermodynamic 
limit — i.e. for small crystallites? Studying the validity 
of the partial differential equations approach at various 
length scales and temperatures is important since this 
formalism is also used by experimentalists to derive dif- 
fusion constants |10 O] or interpret their data |l2]. Some 
workers [^,0 have also used this approach as a black 
box to calculate the time evolution of different structures. 
Our study is related to one of the major problems of equi- 
librium surface physics : the comprehension of the relax- 
ation of a perturbed surface profile, below its roughening 
temperature T^ ||l^. Above T/j, the surface (of an in- 
finite volume crystal whose ratio surface area/volume is 
finite) is rough. This is equivalent to say that the step 
free energy vanishes identically (steps can be created at 
no free energy cost), and it also implies that the surface 
free energy is an analytic function of the local slope. Be- 
low Tr the surface is smooth. This means that the step 
free energy is non- vanishing, and that the surface free 
energy is non-analytic. Indeed, below Tr the surface free 
energy displays cusps at particular orientations, which 
coincide with facets in the equilibrium shape. In the 
thermodynamic limit, the chemical potential of a crystal 
is defined from the Gibbs-Thomson relation, fi — flK, 
in terms of the surface crystal curvature K. On facets, 
the curvature is ill-defined, and the chemical potential is 
fixed by the curved part surrounding the facet. For the 
relaxation of a 2D surface (the surface of a bulk crystal) , 
it is generally believed that MuUins' treatment is correct 
above Tr. Below this temperature, different approaches 
have been proposed [ |l6[ . 

Molecular dynamics (MD) simulations of the coales- 
cence of three-dimensional clusters containing roughly 
1000 atoms, have shown that the relaxation kinetics is 
slower than predicted by MuUins' theory |]1^. However, 
MD simulations are still limited in computation time (no 
more than ~ 10ns) and it is therefore difficult to follow 
the coalescence at temperatures not too close to the melt- 
ing temperature (see also PM]). An alternative method 
consists in using Kinetic Monte Carlo (KMC) simula- 
tions which allow an incomparable larger range of time 
studies. As a first attempt, we have chosen to study a 
two-dimensional (2D) system, namely the relaxation ki- 



netics of 2D faceted islands supported on a triangular 
lattice. 

The basic idea is the following : we start with a island 
with a shape clearly not an equilibrium one (e.g. with 
an X side 10 times longer than the y side) and anneal 
it at a given temperature. Indeed, the perimeter free 
energy dictates the equilibrium shape, which is the one 
that minimizes the island free energy at given volume 
(here we expect hexagonal shapes because of the lattice 
geometry). We assume that the island has relaxed when 
the aspect ratio (defined as the ratio of the x and y gy- 
ration radii, a = Px/Py) becomes 1. We then monitor 
the kinetics of the relaxation process, which depends on 
the precise pathway chosen by the island to change its 
shape. We are particularly interested in the influence of 
facets in the relaxation kinetics. For this reason, we per- 
form simulations at several temperatures: at high tem- 
peratures, where the islands contour is clearly rough and 
only small facets can be distinguished; at low tempera- 
tures, where clear-cut and persistent facets are apparent. 
This point is quite delicate, because a 2D object — whose 
contour is a line, not a surface — is not expected to show 
facets at any temperature in the thermodynamic limit. 
Indeed, a line is rough at all non-zero temperatures, and 
the line tension 7 {Jm~^) is analytic for all orientations 
||15| . However, facets do appear at T = 0, and at low T 
one would expect the persistence length of a facet (the 
average distance between kinks) to be quite large. If it 
is larger than the island side, then facets are indeed ex- 
pected, as we observe. In other words, creating a kink 
costs a finite energy, which is always compensated by the 
entropy gain when the length of the line goes to infinity 
(cf. Landau's argument for the non-existence of phase 
transitions in ID). As long as the line is finite, facets oc- 
cur. One could then guess that they affect the kinetics, 
even for T > Tr, = 0. 

Our main conclusion is that equilibration of an island's 
shape is a non-universal process, in which the time evolu- 
tion of the shape does not obey scaling, while it strongly 
depends on temperature, and thus on system-dependent 
features like the energy scale E. Scaling relations can 
be found for the relaxation (or equilibration) time, as 
a function of temperature and island size. Indeed, two 
regimes with two different scaling forms are born out by 
the simulations, at high and low T, respectively. We ten- 
tatively attribute these two regimes to the absence and 
presence of facets, respectively. 



II. MONTE CARLO AND PARTIAL 
DIFFERENTIAL EQUATION APPROACH 

A. Kinetic Monte Carlo simulations 

We perform "standard" kinetic Monte Carlo (KMC) 
simulations on a triangular lattice. We assume that the 
potential energy of an atom is proportional to its number 



of neighbors, and that the kinetic barrier for diffusion is 
also proportional to the number of initial neighbors, re- 
gardless of the final number of neighbors, i.e. after the 
jump (see Fig. |l|). This is of course a huge simplifica- 
tion, which is however aimed here at describing the global 
evolution of a model island. In other words, we do not 
wish to study any particular system but rather to inves- 
tigate properties which should not depend on the details 
of atom-atom interaction. Therefore, we use a simple 
kinetic model containing as few parameters as possible 
(only one, the ratio E/ksT where E sets the energy scale 
{E — 0.1 eV throughout the paper), ks is the Boltzmann 
constant and T the absolute temperature). Comparing 
with recent ab-initio calculations M for the Al(lll) sur- 
face, we note that our one-barrier assumption does give 
the good order of magnitude of the relative jump frequen- 
cies for the different hopping processes of interest here. 
We also exclude any explicit "Ehrlich-Schwoebel" barrier 
|18| for atoms hopping around corners, although the oc- 
currence of atoms with a single neighbor is treated in a 
special way (see below). The kinetic barriers for some 
jumps are shown in Fig. 0. In the same spirit as ours, 
a similar but slightly more complicated model has been 
used recently by Metiu's group with the scope of obtain- 
ing system-independent information on island diffusion 
on a surface: these authors investigate the existence of a 
"universal" size dependence of the island diffusion con- 
stant. They conclude that such universality is not ob- 
served, and we observe a similar phenomenon for island 
equilibration. 

The time evolution of the island shape is obtained 
by the following algorithm. We first calculate the 
following quantities : Si — exp[— i?/(10fcBT)], S2 — 
exp[-2E / (kBT)] and S3 = cxp[-3E / (kBT)] which rep- 
resent the relative weights for the jump probabilities for 
atoms with respectively 1, 2 or 3 neighbors (atoms with 
more neighbors simply do not move : see an explanation 
of the precise forms of the different Si below). Then, in 
each iteration, we calculate the probability to move an 
atom with i neighbors as : 



(6 - i)ni5i 



(2.1) 



where Ui is the total number of atoms having i neigh- 
bors. We choose randomly one of the atoms with the 
appropriate number of neighbors and move it in a ran- 
dom direction. The time is increased at each iteration 

by 



dt 



3 



i^o > (6 - i)niSi 



I / 



(2.2) 



where vq is a Debye frequency (we have taken vq = 
10^^ s~^). To check that the law of detailed balance is sat- 
isfied, one can refer to Fig. |l| : the probability for an atom 
to jump from a site having n neighbors to a site having p 



neighbors is 6n while the opposite transition has a proba- 
biUty 6p. Their ratio is 5n/Sn = exp{—E/{kBT){n — p)), 
i.e. equal to the energy difference of the initial and final 
configurations, as required by the law of detailed balance 
(the particular case when n or p is equal to 1 can be 
analyzed in the same way). This algorithm is very fast 
[jl9|-pl[ since all iterations contribute to the evolution 
(there are no rejected moves). One peculiarity of this 
model is the treatment of atoms having one single neigh- 
bor: (5i is much larger than what one could expect from 
the general rule (5„ = cxp[—nE/{kBT)]: indeed, we let 
(5i = exp[-^/(10fcBT)] instead of Ji = exp[-^/(fcBr)]. 
This is to ensure that singly-bonded atoms, which are in 
some sense in a "transition state" , rapidly go into some 
physically reasonable position, i.e. one having 2 or more 
neighbors. Note also that detachment of atoms from the 
islands is forbidden here: equilibration is only due to 
mass transport along the island contour. This is clearly 
different from Ostwald ripening where islands evolve in 
equilibrium with a two-dimensional adatom gas. We note 
that a recent experimental study by Stoldt et al. |g^ has 
shown that supported Ag two dimensional islands do in- 
deed relax via atomic diffusion on the island edge, with- 
out significant contribution from exchange with a two- 
dimensional adatom gas. A last remark on the algorithm 
used here : we do not allow atoms having more than 
3 neighbors to move. In some sense, they have an infi- 
nite potential energy. Since our potential energy is not 
very realistic anyway, this hypothesis allows to simplify 
and accelerate the simulations. The key point is that our 
results are particularly interesting at low temperatures, 
where including the possibility for atoms with 4 neigh- 
bors to jump would make no significant difference in the 
kinetic evolution of the island because at these tempera- 
tures their jumping is vanishingly small. 



B. Partial Differential Equation approach 

A complementary approach at predicting the evolu- 
tion of a crystal shape at a temperature higher than 
the roughening transition, consists in coarse-graining the 
crystal profile, in order to treat it as a smooth function 
/i(x, t), and in writing down its time evolution in the form 
of a partial differential equation, whose form depends on 
the physical situation of interest. The situation when 
matter transport is assured by adatom diffusion along 
the surface, has been originally considered by MuUins 
and coworkers [p[ , for studying small deformations of an 
infinite planar surface. The case of a finite, closed "sur- 
face" — the island contour — is somewhat more subtle, 
and we give the derivation in some detail. The evolution 
equation has in general then the form of an equation for 
s the curvilinear coordinate or arclcngth. On purely geo- 
metrical grounds, it can be shown p3| that this equation 
can be written as an evolution equation for the curvature 
K{s, t), of the form 



dK 






K' 



dK 
ds 



ds'Kvn (2.3) 



where Vn is the normal velocity of the interface. The 
latter is fixed by the physics of the problem. In our case, 
when edge diffusion is the relevant physical process de- 
termining the relaxation of the shape, the equation for 
Vn must have the form of a conservation equation for the 
island area 



3. 

ds 



(2.4) 



The edge diffusion current j is given by the gradient of 
the local excess chemical potential as in Pick's law : 



J = 



D dAn 
ksT ds 



(2.5) 



where D is a (collective) perimeter diffusion coefficient. 

The excess chemical potential is in turn related to the 
local curvature K{s) through the Gibbs-Thomson rela- 
tion: 



A^ = 



-ja^K 



(2.6) 



where 7 is the line tension (that for simplicity is as- 
sumed isotropic here), and a the depth of the outer layer 
within which mass transport takes place (one lattice spac- 
ing in our case). Eqs. (p.Sj) and (2.6) yield 



DjadK 
knT ds 



(2.7) 



and the corresponding evolution equation (2.4) 
given, in our two-dimensional situation, by: 



IS 



Dja^ d^K 
ksT ds^ 



(2., 



Equations (2.3) and (2 



must be solved simultane- 
ously f or the island shape: the results will be shown in 
section HID . The equilibrium solution w„ = obviously 
is a constant-cu rvat ure shape^ that is a circle. Also, note 
that equations (2.3) and ( p.8|) are invariant by a rescal- 
ing s — > As, i — > X*t M, so that the equilibration time 
of a deformed island of size L is expected to be propor- 
tional to L^ or N^ where N is the number of atoms in 
the island. 



III. MONTE CARLO SIMULATION RESULTS 

To study the influence of the facets on the coalescence 
kinetics, it is interesting to study the time evolution of 
the aspect ratio and the size, as well as the temperature 
dependence of the equilibration time teq- We recall that 
teq is defined as the time needed for the island to reach 
its equilibrium shape. In practice, we take teq as the first 
time when the aspect ratio a defined above becomes less 
than 1. Each point is the average of several runs (up to 
200 for the smallest islands). 



A. Island morphology 

Fig. m shows the time evolution of the perimeter of 
a 6250 atoms island at 500 K {E/ksT = 2.3). It is 
clear that the shape evolution occurs with rough island 
borders. 

Fig. yb shows the time evolution of the perimeter of 
an island containing 6250 atoms at 83 K [E/ksT = 14). 
At this temperature facets are apparent throughout the 
evolution. 

A more precise comparison of the presence of facets at 
the two temperatures studied above is given in Fig. ^ It 
is apparent that facets are present at 83 K, in contrast 
to the high rugosity observed at 500K. 



B. Dependence of the equilibration time on island 
size 



3 and 4 times the energy needed to break a single bond. 
A tentative interpretation of these values is given below 
in Sec. Irvl. 



The continuous analysis of section II B predicts teq ^ 
N"^ as a function of the number of atoms N inside the 
island, for any temperature. Indeed, the numerical so- 
lution of the full, non-linear equations ( |2.3| ) and ( |2.§| ) 
appears to agree with this prediction. Fig. 0a shows 
the size dependence of t^q for different temperatures as 
given by the simulations. The simulation results agree 



with te 



N only at high temperatures. Below 250 



K, it is clear that teq increases slower than A^^, and the 
lower the temperature, the smaller the exponent. One 
can also notice (Fig. ||b) that the local exponent for low 
temperatures approaches 1 for the highest island sizes. 
This is analyzed in Sec. |V| where we give an attempt 
at deriving a scaling relation describing the two regimes. 
It should also be noted that extrapolating the different 
curves for very high values of the island size leads to 
an apparently absurd conclusion : very large islands do 
equilibrate faster at lower temperatures. This is a imme- 
diate consequence of the higher size exponents found for 
the highest temperatures. To avoid a paradox, we must 
admit that there exist a crossover from high to low tem- 
perature behavior for a given size that depends on the 
temperature. Therefore, even at 83 K (highest curve), 
for large enough islands, one should recover the t^q ~ A^ 
regime. The scaling analysis presented below explains 
this crossover. 



C. Dependence of the equilibration time on 
temperature 

Fig. Ba shows that t^q rapidly increases as tempera- 
ture decreases, in roughly the same way for all the island 
sizes. The equilibration time is not exactly a thermally 
activated quantity, since there is a clear curvature in its 
Arrhenius plots, as shown in Fig. gb: the local activa- 
tion energies increase from roughly 0.3 eV at high T to 
0.4 eV at low temperature. This represents respectively 



D. Precise kinetics of the relaxation 

We have seen that the size dependence of the equili- 
bration time obeys the scaling predicted by the linearized 
equations only at high temperatures. It is interesting to 
chec k whether the full solution of equations ( p.3| ) and 
( ^.8[ ) agrees with the high-temperature behavior of the 
MC simulations, when no facets are apparent, and the 
island looks rough and rather isotropic. Fig. |^, where 
the aspect ratio is plotted as a function of the reduced 
time t/teq, shows that no agreement is found, at any of 
the studied temperatures. Indeed, this is a posteriori not 
surprising, since the MC results do not seem to obey any 
scaling relation, or maybe only at high temperature, and 
it is then obvious that the "universal" description given 
by continuous equations does not apply. It is nevertheless 
a little surprising that the continuous description, which 
agrees with simulations in the case of a planar surface 
above the roughening temperature, does not seem to set 
a limiting behavior valid for very large sizes (A — > cxd), 
nor to provide the scaling form which seems to appear in 
the simulation results at high temperatures. 

We can think of (at least) three explanations of this ob- 
servation. First, it could be argued that this is an effect 
of the edge tension. Indeed, in writing the constitutive 
equation (2.8) we have assumed that 7 is isotropic. This 
is clearly not the case in the simulations since we take a 
triangular lattice, and the energy of the facets depends 
on the orientation. We are currently performing a nu- 
merical integration of equations (p]q) and ( ^.8[ ) including 
anisotropic line tension to further investigate this point. 

Second, one could question the adequacy of the con- 
tinuous treatment to describe the detailed path to equi- 
librium, even at high temperatures, for small clusters. 
Indeed, one could argue that the macroscopic concepts 
of curvature, chemical potential, etc. are not adapted to 
deal with nanometric objects containing few (less than 
10000) particles. 

Finally, it is possible that the simulations are not 
adapted to agree with the continuous theory, in the sense 
that our interatomic potential is too crude to give a rea- 
sonable kinetic path to equilibrium. It is clear that the 
assumption that the transition probability depends only 
on the initial state is n ot ge nerally correct. As has been 
argued above (section II A ), one would expect such a 
rough potential to reproduce a universal exponent (as 
is observed for tf,q oc A^ at high temperatures) but not 
necessarily a detailed time evolution, if the latter is non- 
universal. 



IV. A SCALING ARGUMENT 

Our MC results show that two different regimes - at 
low and high temperatures - can be identified. At first 
sight, this is surprising since the only occurrence that 
could separate high from low temperatures is the rough- 
ening transition, which, strictly speaking, takes place at 
T =0 K in 2 dimensions — that is, for a one-dimensional 
"surface". However, facets do not disappear suddenly as 
the temperature is raised from T = 0. Indeed, a "persis- 
tence length" of facets can be defined as the equilibrium 
value of the distance between kinks along the step edge. 
At equilibrium and at low temperatures, we can consider 
an "ideal gas" of kinks, whose density (number of kinks 
per unit length) is given by the formula 124] 



rikink = 2eyi^{-(3W)/a 



(4.1) 



where W is the kink creation energy and a the lattice 
spacing. Indeed, to form kinks one has to take an atom 
out of the step edge, and to place it anywhere along the 
step. On a triangular lattice, in doing so one looses 4 
nearest- neighbours bonds, and gains 2. A net balance 
of 2 broken bonds results. In the process, 4 kinks have 
been created (each atom counts for 2 kinks), so that the 
energy cost per kink is 



W = E/2. 



(4.2) 



The factor of 2 in ( |4.l| ) comes from the fact that kinks 
always appear and disappear in pairs — in other words 
two types of kinks, positive and negative exist, of equal 
number. Then, the equilibrium distance between kinks 
is 



£o = aexp(/3£;/2)/2, 



(4.3) 



which diverges as T goes to 0. At a given temperature, 
a step looks straight (free from kinks) over lengths of the 
order of Iq. An island is thus bound to look faceted as 
long as ^0 is larger than the island linear size L. Then, 
the approximate equality 



exp(£;/(2fcBT,)/2 « Lja 



(4.4) 



gives the "crossover size" L^ (at fixed temperature) 
or the "crossover temperature" T^ (at fixed size) for the 
crossover between the high (rough) and low (faceted) 
temperature regimes. A comparison with Fig. (5b), where 
the high temperature regime corresponds to teq ~ N"^, 
while the low temperature regime corresponds to tf,q ^ 
N, shows that this criterion is not too bad: the formula 
predicts a crossover temperature for an island of size 
N — 500 {L w 22) of approximately 240 K, in good agree- 
ment with the simulation data. Note that for smaller 
sizes than about N « 100 the simulations do not show 
a well defined low temperature behaviour. We attribute 
this to the large importance of geometric kinks (imposed 
by the fact that the step closes on itself) over thermal 



kinks for these small sizes. Indeed, we claim that the low 
temperature regime is ruled by the equilibration of the 
spatial distribution of thermal kinks: the initial shape 
creates a strongly inhomogeneous distribution of these 
kinks, which then diffuse to achieve spatial uniformity, 
and thus equilibrium. Kinks diffuse by emitting atoms, 
so that we conclude that atom emission from kinks is the 
limiting kinetic step determining the low temperature be- 
haviour of teq ■ 

Based on this assumption, we can give a scaling argu- 
ment that reproduces well the observed t^q as a function 
of N and T. The argument is similar to that used by 
Bales and Zangwill p^] and Pimpinelli et al. p^ to dis- 
cuss step roughening and smoothening during growth and 
at equilibrium. Indeed, it amounts to performing a linear 
stability analysis, and computing the relaxation time of 
a perturbation of given amplitude and wavelength. For 
sake of clarity, we start from the discussion of the high 
temperature behaviour, which is well known from MuUins 
work. 

Physically, we assume that there are two ingredi- 
ents determining the relaxation : the thermodynamical 
"force" which drives the relaxation (here, the excess cur- 
vature) and second, the kinetic factors which determine 
the rate of the equilibration. At all temperatures, curva- 
ture effects are relevant, but we assume that the kinetics 
change due to the presence of facets (or, equivalently, 
to the low concentration of kinks). The transition takes 
place, as stated above, when io ~ Lc- Here is the math- 
ematical translation of this idea. 



A. High temperatures 

Let 5q be the amplitude of a perturbation of wave vec- 
tor q of the island perimeter with respect to the equi- 
librium shape. The curvature effect (Gibbs-Thomson) 
opposes the increase of the deformation. The rate of de- 
crease depends on the appropriate kinetic process which 
limits transport of matter from high to low chemical po- 
tential regions. Let n^q be the equilibrium atom density 
along a reference island edge with the equilibrium shape. 
Then, at high temperature, a deformation of local cur- 
vature K results in an e xcess chemical potential A// ^ 
^q^5q as in equation (^^). In turn, this creates an excess 
atom density rieajc = Ue, exp[A^/(fcBT)] « neq[l+Tq^5q], 
where F = ^/{kBT). Then, edge atoms flow away from 
the deformation, whose amplitude decreases at a rate 
proportional to the divergence of the mass current : 

'^1 ~ -^V2(ne,c - neq) « -^92 ^ pg2^^ (4 5) 

where r* is the typical timescale of the appropriate ki- 
netic process which is responsible of matter transport. 

A more detailed justification of this expression can 
be found in MuUins [], Bales and Zangwill |2^ and 
Pimpinelli et al. p6|]. 



Defining the equilibration time teq by writing (5q=i/L 

-Sq^l/L/teq giveS 



4- « Dnl = ,y^a/DoeM~AE/{kBT)] 



(4.11) 



i-eq 



«egr 



iV^ 



"e^r 



(4.6) 



Inserting ( [4.1l|) and (^ 
temperatures, 



in Eq. 4.10 yields, at low 



MuUins equation is recovered if one assumes that atom 
edge diffusion limits the kinetics, so that 



teq^^Nexp[4E/{kBT)] 



(4.12) 



Again, the activation energy predicted here is in good 
D. (4-7) agreement with the low temperature limit observed in the 

simulations (Fig. g). The scaling teq ~ iV is less clearly 

seen in the simulations (Figs. |^). However, the simula- 

where tions show that the lower the temperature, the lower the 

i/Q exp [-3E/{kbT)] is the rate of atom emission size exponent, and if A^ is not too small, teq - iV is con- 



1 



The atom equilibrium density can be obtained from the 



^kink , 



detailed balance at the kinks: Due 

^kink 

from kinks and Dricq is the atom flux to the kinks 
Thus, 



''eq 



= iyo/D„exp[-E/{kbT)]. 



(4- 



Inserting (17) and (| 
of high temperatures. 



^) in Eq. |4.6| yields, in the limit 



eq 



1 



N'^e^p[3E/{kbT)] 



This prediction reproduces the t, 



eq 



(4.9) 



N'^ scaling of 



the continuum theory, and it is in very good agreement 
with the simulation results obtained at high tempera- 
tures both for the temperature dependence and for the 
size dependence (Fig. ^. Indeed, at high temperature 
the equilibration time shows an activation energy of ap- 
proximately 3E (Figs. H), 
regime. 



and teq behaves 



N^ 



this 



B. Low^ temperatures 

The low temperature regime sets in, for a given crystal 
size, when the equilibrium distance between kinks be- 
comes of the order of the linear size of the crystal, and 
straight step portions appear. The (thermal) kink den- 
sity then becomes a relevant concept. When the crystal 
is deformed from the equilibrium shape, the kink density 
is increased where the facets are shrunk, and decreased 
where they are streched. On removing the constraint, 
the kink density tends to equilibrium and seeks spatial 
uniformity. If the equilibrium facet size is L « A^^/^, 
and a shape deformation of order Si is introduced, the 
kink density unbalance is approximately SijL?. Then, 
the perturbation relaxes as 

{k) « — ^ X \bl « ^-\-M. (4.10) 

The relaxation proceeds by moving a whole row of atoms 
from a short to a long facet; diffusion is fast on facets, 
and the process is limited by nucleation of the new row, 
that is, by the rate of atom encounters Dn^ . Then, 



sistent with our results. When N is smaller than about 
100, teq seems to increase faster than linearly. At such 
small sizes, facets are always very short, and it is likely 
that an intermediate behaviour between mass transport 
and facet nucleation rules the relaxation. 



C. Discussion 

The scaling argument we propose nicely reproduces 
the results of our simulations and leads to a reasonable 
physical picture of the equilibration, consistent with the 
observed morphologies and kinetics (presence of facets, 
rapid completion of atomic rows . . . ). 

Even more, our results can be used to estimate the 
behaviour of the diffusion coefficient of a cluster as a 
function of the cluster size N k. 1? ,\>y means of another 
scaling argument. In order to diffuse over a length £, 
a number of atoms of order IL have to be transferred 
from one side of the island to the opposite side. The 
time needed to do this is of the order of the time teq{L) 
needed to equilibrate a fluctuation of linear size L and 
mean square amplitude £^ « L/{P"/) ||l5[. Therefore, 
from the knowledge of teq{L) we can know the diffusion 
coefficient D{N) from the Einstein relation 



D{N)teq{N) ^f -7Vi/2_ 



(4.13) 



If we assume that our high-temperature result 
teq{N) - 7V2 holds, we find 



D{N) - N-^/'^. 



(4.14) 



If we assume that our low-temperature result teq{N) 
N holds, we find 



D{N) - Ar-i/2. 



(4.15) 



Equations (4.14) and (4.15) can be compared with the 
results of the simulations of island diffusion of Metiu 
and coworkers': on a (OOl)-type substrate they find that 
the size-dependent diffusion constant D{N) of 2D islands 
varies as D{N) ~ 1/iV^-^^ at high temperature, and as 
D{N) ^ l/7V°-^2 at low temperature M. Of course. 



different equilibration processes would lead to different 
teq{N)s and then to different behaviours for D{N). This 
might explain the different results obtained by Bogicevic 
et al. |3y] for islands diffusing on a (111) substrate but 
with energy barriers for the jumps different from those 
assumed here. 



V. SUMMARY, PERSPECTIVES 

The relaxation to equilibrium of 2D islands contain- 
ing up to 20000 atoms shows unexpected features. Our 
results show that there is no "universal" size exponent 
for island equilibration, a result similar to that found by 
Metiu's group for island diffusion [^. We are now study- 
ing the case of 3D clusters to check both the scaling of 
the equilibration time with the size of the particle and the 
precise kinetic path followed to reach equilibrium. This 
is done by KMC simulations and an analytical approach. 



e-mail address : jensen@dpm.univ-lyonl.fr 
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(a) 



FIG. 1. Potential energy of an atom diffusing along the is- 
land edge when it has 1 (energy : -E/10), 2 (energy : -2E) or 
3 (energy : -3E) neighbors. The atom energy is supposed to 
depend only on the number of first neighbors and the transi- 
tion state is assumed to lie at the same energy for all jumps, 
taken here as the origin of energies. As a consequence, the 
energy barriers for diffusion are equal to the potential energies 
of the atoms, ensuring detailed balance (see text). 
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(b) 
FIG. 3. Time evolution of islands containing 6250 atoms 

at 500 K (a) and 83 K (b). The initial state corresponds to 

the most elongated configuration. 
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(b) 
FIG. 2. Examples of activation energies for different atomic 

jumps on the island edge. Note that detachment of atoms 

from the island is explicitly forbidden in the simulations. 




FIG. 4. Large facets are present at 83 K (open circles) in 
contrast with the rugosity observed at 500 K (filled triangles) 
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FIG. 6. (a) Time needed to reach equilibrium as a function 
of the temperature for different island sizes, (b) shows the 
local activation energy (defined as in Fig. H) at each temper- 
ature. 



FIG. 5. (a) Size dependence of teq for different tempera- 
tures and (b) 
local slope of the size dependence (obtained as the discrete 
derivative of the ieg vs A'' curve given in (a) : for point i, it is 
{log{t,qii+l))-logiUqii-l)))/{log{N{i+l))-log{N{i-l))). 
The local slope remains close to 2 at high temperatures for all 
sizes, but it approaches 1 for high sizes at low temperatures. 
In (a), the curves have been shifted vertically for clarity. The 
precise fits are the following : T=500 K : teq = 8 W~'^^N^-°°; 
T=250 K : teq ^ 8 10'^ N^-^^; T=125 K : t^q = CSIiV^'^^- 



- Mulllns 

• 6250 T 500 
6250 T 250 

n 6250 T 167 
T 6250 T 100 

♦ 6250 T 83 

- 640 T 83 




T=100 K:teq = 6 IQ^N^-^^; T=83 K : te 



= 9 10'' iV^-^^ 



FIG. 7. Kinetic path followed by the islands as they re- 
lax to equilibrium. The solid line represents the prediction of 
MuUins' equation 2.i and the different symbols to the relax- 
ations obtained in the KMC simulations for islands containing 
6250 atoms at the temperatures shown. The dashed line cor- 
responds to a smaller island and illustrates the wide range of 
relaxation paths observed. 
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